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Department  of  Metallurgical  and  Materials  Engineering 
University  of  Pittsburgh 


ABSTRACT 

This  report  summarizes  work  on  composite  and  alloy  solidification 
presently  conducted  under  the  ONR  contract  and  includes  two  papers  on  the 
simulation  of  heat  flow  and  thermal  stresses  in  continuous  casting 
presented  at  the  "Conference  on  Computer  Simulation  for  Materials  Applications" 
in  Gaithersburg,  Md.  in  April  1976. 

A finite  element  analysis  of  heat  transfer  in  continuous  casting  is 
presented  in  the  first  paper.  The  steady  state  conditions  of  the  solidifica- 
tion of  a cylindrical  alloy  ingot  are  simulated.  The  analysis  takes  into 
account  both  radial  and  axial  heat  conduction  and  the  variation  of  alloy 
properties  with  temperature.  The  analysis  has  the  versatility  to  simulate 
any  casting  configuration  as  the  program  has  the  capability  to  treat  any 
arbitrary  boundary  conditions.  In  the  mold  region,  where  an  air  gap  generally 
forms,  the  position  of  gap  formation,  the  size  of  the  gap  and  the  heat 
transfer  through  the  gap  by  conduction  and  radiation  are  computed.  Values 
predicted  bv  the  model  are  in  close  agreement  with  published  Laboratory 
and  industrial  data. 

The  second  paper  presents  a three  dimensional  analysis  of  steady  state 
thermal  stresses  developed  in  the  continuous  casting  of  cylindrical  sections. 

The  analysis  takes  into  account  the  elastic,  plastic  and  steady  state  creep 
deformations  of  the  castings.  It  can  take  into  account  virtually  all 
boundary  conditions  encountered  in  continuous  casting  (friction  in  the  mold, 
metallostatic  head,  etc.)  and  is  applicable  to  both  ferrous  and  nonferrous 
continuous  casting.  The  results  from  the  analysis  are  compared  to  the 
published  residual  stress  measurements  made  on  aluminum- 12%  silicon  alloy  (A132I. 
The  model  is  very  useful  in  understanding  the  influence  of  the  operating 
parameters  on  the  hot  tearing  and  cold  cracking  of  continuous  castings. 
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SUMMARY  OF  WORK 


Ji 

The  research  supported  under  this  ONR  contract  has  concentrated  on  two 

overlying  technological  themes: 

Control  of  structure,  microsegregation,  and  macrosegre- 
gation in  ingots  and  castings  by  control  of  solidification 
parameters . 

Development  of  high  temperature,  high  strength  composite 
materials  by  solidification  techniques. 

Work  on  melt  grown  conposite  materials  has  been  with  eutectic  and 
peritect ic  alloys.  A continuing  theme  in  eutectics  has  been  establishing 
the  relation  between  growth  conditions  in  zone  melting  and  microstructure  for 
Pb-Sn,  a low  temperature  model  system,  and  Nb-Nb2C,  a high  temperature  model 
system.  Present  work  in  this  area  deals  with  defining  changes  in  the  micro- 
structure that  occur  immediately  behind  the  growth  interface  in  Nb-Nb^  and 
determining  whether  there  is  a plate  to  rod  change  in  morphology  with  growth 
rate  in  Pb-Sn.  Such  a transition  in  morphology  with  growth  rate  has  been 
demonstrated  for  Nb-Nb2C. 

Studies  on  the  high  temperature  properties  of  Nb-Nb2C  have  been  undetaken 
with  the  aim  of  defining  the  effect  of  microstructural  variables  on  the 
properties,  thermal  stability  and  creep  resistance.  Nb-Nb2C  composite  samples 
with  plate  and  rod-like  carbide  morphologies  showed  good  microstructural 
stability  when  exposed  to  temperatures  up  to  1125°C.  On  exposure  of  these 
samples  to  temperatures  in  excess  of  one  half  the  eutectic  temperature,  the 
structures  coarsen  appreciably.  For  samples  with  plate  carbide  morphology, 
the  plate  morphology  and  high  aspect  ratio  is  maintained.  For  samples  with 
rod-like  carbide  morphology,  the  major  change  was  a reduction  in  the  total 
number  of  carbide  fibers  and  a transition  in  carbide  morpholo©'  from  rods  to 
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plates.  Present  optical  and  S.E.M.  observations  indicate  the  change  in  morphology 

on  annealing  is  a two  dimensional  coarsening  with  preferential  growth  of 

the  surviving  carbide  rods  on  the  (100)  faces  until  they  merge  to  form  plates. 

Emphasis  in  the  creep  work  has  been  placed  on  relating  the  microstructue 
of  eutectic  Nb-Nb2C  composites  to  their  creep  properties.  It  is  worth  noting 
that  the  steady  state  creep  rate  at  1100°C  of  well  aligned  Nb-Nb2C  composites 
is  the  order  of  the  creep  rate  of  Ni-base  superalloys  at  800°C.  There  is  a 
very  strong  dependence  of  creep  rate  on  structure.  The  creep  rate  can  increase 
by  two  orders  of  magnitude  in  going  from  a well  aligned  composite  structure 
to  a well  aligned  dendritic  structure.  Misorientation  and  banding  have  a 
major  influence  on  creep  rate.  The  stress  exponent  varies  from  3 to  14 
depending  on  microstructure.  The  dependence  of  activation  energy  on 
microstructure  (if  any)  is  being  measured.  The  metallographic  investigation 
of  samples  after  creep  testing  is  underway. 

The  first  work  on  peritectic  alloys,  the  Pb~Bi  system,  brought  surprising 
results.  Aligned  two  phase  composite  structures  were  produced  and  these  were 
assumed  to  be  the  result  of  dendritic  growth  at  high  degrees  of  undercooling. 

It  was  also  possible  to  produce  single  phase  supersaturated  solutions  in 
compositions  that,  at  equilibrium,  should  be  a two  phase,  peritectic  mixture. 

We  have  now  confirmed  the  results  for  Sn-Cd.  New  work  will  be  on  higher 
temperature  peritectic  systems. 

The  work  on  producing  an  oxide  dispersion  strengthened  alloy  by  casting 
has  met  with  mixed  results.  We  have  been  able  to  rheocast  vacuum  melted 
alloys  and  to  mix  in  very  fine  oxide  dispersions.  However,  we  have  not  been 
able  to  mix  and  cast  sufficiently  quickly  to  avoid  reaction  between  the  oxide 
and  solute  in  the  melt.  Proposed  work  would  be  small  scale  experiments  to 
measure  the  effect  of  surface  treatment  and  other  particle  variables  in 
delaying  the  reaction. 


3. 

A computer  program  has  been  developed  to  simulate  the  steady  state 
continuous  casting  of  round  ingots.  The  program  calculates  temperature 
distributions  throughout  the  ingot  and  build  up  of  thermal  stresses  as  a 
function  of  casting  variables  and  alloy  properties.  Presently,  parametric 
studies  are  being  simulated  to  see  the  individual  and  combined  effects  of 
the  variables.  The  next  sections  of  this  report  deal  with  the  analysis 
method  of  the  continuous  casting  simulaton. 

The  remaining  alloy  work  deals  with  titanium  and  N1TIN0L.  The  titanium 
work  will  evaluate  the  tendency  of  these  alloys  to  exhibit  macrosegregation 
as  a result  of  the  bulk  movement  of  enriched  liquid.  The  NITINOL  work  will  first 
aim  at  developing  a predictive  ability  with  regard  to  manufacturing  process 
and  transition  range  and  second  will  be  to  prepare  a modest  inventory  of 
research  materials. 

COMPUTER  SIMULATION  OF  CONTINUOUS  CASTING 

Much  of  the  basic  work  we  have  carried  out  under  this  contract  is  applicable 
to  directional  solidification  processes;  that  is,  casting  methods  wherein  there 
is  a positive  temperature  gradient  in  the  liquid  ahead  of  the  dendrite  tips. 
Continuous  casting,  electroslag  remelting,  and  vacuum  arc  remelting  are 
examples  of  directional  solidification  processes.  A comprehensive  engineering 
model  of  these  processes  would  include  being  able  to  predict  heat  flow  patterns, 
fluid  flow  patterns,  states  of  stress,  grain  and  dendrite  structure,  micro- 
and  macrosegregation,  inclusion  shapes  and  distributions,  hot  carcking,  and 
residual  stresses.  One  aspect  for  which  we  have  developed  a model  is  the 
build  up  of  thermal  stresses  in  continuous  casting. 

The  numerical  simulation  of  continuous  casting  has  been  developed  by 
J.  Mathew  for  his  Ph.D.  thesis.  The  simulation  utilizes  the  finite  element 
method  and  considers  the  steady  state  portion  of  the  casting  of  round  ingots. 

In  order  to  make  the  thermal  stress  analysis  of  continuous  casting,  it  is 
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prerequisite  that  analysis  be  made  of  the  steady  state  temperature  distribution 
in  the  various  regions  of  the  ingot  and  certain  properties  of  the  alloy 
being  studied  be  available  as  functions  of  temperature.  The  important  parameters 
include  phase  relations,  stiffness,  relaxation  characteristics,  thermal 
expansion  coefficient,  yield  strength,  and  solidification  characteristics. 

In  order  to  make  the  simulation  appropriate  nonferrous  as  well  as  ferrous 
casting,  both  radial  and  axial  heat  transfer  are  taken  into  account,  i.e., 
the  heat  flow  analysis  is  two  dimensional.  The  stress  analysis  is  three 
dimensional  and  has  the  capability  of  accounting  for  elastic,  plastic,  and 
creep  strains.  The  inclusion  of  creep  behavior  in  the  simulation  reflects  our 
feeling  that  the  build  up  of  thermal  stresses  is  strongly  dependent  on  time  and 
rate.  The  details  of  the  analysis  are  covered  in  the  papers  included  with  this 
report . 

Presently,  the  program  is  being  applied  to  parametric  studies  of  the 
effect  of  casting  variables  and  alloy  variables  on  the  thermal  pattern  and 
the  build  up  of  thermal  stresses  in  continuotfi  casting.  The  casting  parameters 
being  considered  are  section  size,  mold  length,  cooling  sequence  and  intensities, 
casting  rate,  casting  superheat,  and  pressure  head.  The  program  is  also 
being  applied  to  several  actual  cases  where  alloy  and  casting  data  are 
available.  As  is  covered  in  the  papers,  the  results  of  the  program  are  in 
good  agreement  with  available  published  data.  Another  aspect  of  the  present 
effort  is  the  documentation  of  the  program  to  the  state  it  can  be  widely 
used. 

The  present  program  has  been  built  in  modular  fashion  so  that  in  subsequent 
phases  of  the  development  we  will  be  able  to  use  the  existing  program  as  a 
framework  for  simulating  the  vertical  casting  of  slabs  and  the  influence  of 
perturbations  in  casting  conditions. 


Presented  at  Conference  on  Computer  Simulation  for  Materials  Applications, 
Gaithersburg,  MD,  April  1976. 

ANALYSIS  OF  HEAT  TRANSFER  IN 
CONTINUOUS  CASTING  USING  A 
FINITE  ELEMENT  METHOD 

J.  Jfcthew  and  H.  D.  Brody 
Metallurgical  and  Mterlals  Engineering 
University  of  Pittsburgh 


Abstract 

A finite  element  analysis  of  heat  transfer  in  continuous  casting  is  presented.  The 
steady  state  conditions  of  the  solidification  of  a cylindrical  alloy  ingot  are  simulated. 
The  analysis  takes  into  account  both  radial  and  axial  heat  conduction  and  the  variation 
of  alloy  properties  with  temperature.  The  analysis  has  the  versatility  to  simulate  any 
casting  configuration  as  the  program  has  the  capability  to  treat  any  arbitrary  boundary 
conditions.  In  the  mold  region,  where  an  air  gap  generally  forms,  the  position  of  gap 
formation,  the  size  of  the  gap  and  the  heat  transfer  through  the  gap  by  conduction  and 
radiation  are  computed.  Values  predicted  by  the  model  are  in  close  agreement  with 
published  laboratory  and  industrial  data. 

Introduction 

This  paper  presents  part  of  a simulation  program  (PCCTSS-Pitt  Simulation)  developed  at 
the  University  of  Pittsburgh  to  study  the  continuous  casting  process.  This  paper 
describes  the  heat  transfer  analysis.  An  accompanying  paper  (1)  describes  how  the  data 
from  the  heat  transfer  analysis  can  be  used  to  compute  the  three  dimensional  stresses 
and  strains  in  the  casting.  The  two  sioulations  are  Intended  to  provide  a valuable  tool 
for  the  study  and  design  of  cylindrical  continuous  casting  processes. 

The  present  thermal  analysis  is  developed  so  that  It  can  simulate  both  ferrous  and 
nonferrous  continuous  casting.  In  Aluminum  D.C.  casting,  the  high  conductivity  of  the 
material  and  the  relatively  slow  cooling  speeds  employed  make  the  axial  heat  transfer 
relatively  important.  In  ferrous  continuous  casting,  axial  heat  transfer  may  be 
neglected  except  for  the  mold  exit  region  where  the  Increased  radial  resistance  to  heat 
flow  caused  by  the  at r gap  formation  make  the  axial  heat  transfer  relatively  important. 
None  of  the  published  analyses  of  heat  transfer  in  continuous  casting  to  this  date  (2-15, 
17-21)  (except  for  reference  IS)  takes  into  account  the  axial  heat  transfer  in  the  cast- 
ing. The  analysis  presented  here  accounts  for  the  axial  and  radial  heat  transfer  in  the 
cylindrical  section,  the  variation  of  thermophysical  properties  of  the  material  with 
temperature,  and  the  air  gap  formation  in  the  casting.  In  the  air  gap,  heat  transfer 


du«  to  conduction  and  radiation  art  considered.  The  coupling  effect  of  the  air  gap 
formed  and  the  heat  transfer  are  simulated  by  resorting  to  an  iterative  scheme.  The 
versatility  of  the  finite  element  method  used  in  this  analysis  enables  the  use  of  any 
arbitrary  boundary  condition  anywhere  in  the  casting. 

The  previous  investigations  of  heat  transfer  in  continuous  casting  may  be  described  under 
(a)  analytical  and  empirical  methods  (3-7)  (b)  integral  profile  method  (8-15)  and  (c) 
numerical  (finite  difference)  methods  (13,16-21).  The  empirical  methods  are  very  limited 
in  scope  of  application.  The  integral  profile  methods  are  very  simple  to  use  and  give 
fairly  good  results  when  applied  to  the  continuous  casting  of  steel.  However,  the 
integral  profile  methods  do  not  account  for  the  heat  transfer  in  the  axial  direction  in 
the  casting  or  for  the  variation  of  thermo physical  properties  of  the  material  with 
temperature.  The  finite  difference  methods  are  very  powerful  and  have  been  shown  to 
give  reliable  results  when  applied  to  the  continuous  casting  of  ferrous  and  nonf»r-...us 
metals. 

The  present  analysis  employs  the  finite  element  method  for  the  simulation  of  heat  trans- 
fer in  continuous  casting.  In  this  method,  the  casting  is  represented  by  a system  of 
triangular  ring  elements.  An  approximate  solution  for  the  temperature  field  is  assumed 
within  each  element  and  the  heat  flux  equilibrium  equations  are  developed  at  a discrete 
number  of  points  within  the  elements.  These  equations  are  assembled  and  solved  for  the 
temperature  distribution  within  the  casting.  This  method  is  general  and  is  shown  to  be 
applicable  to  a wide  variety  of  problems  (22).  The  finite  element  heat  transfer  analysis 
presented  here  has  all  the  advantages  of  the  implicit  finite  difference  methods.  The 
program  is  completely  automated  and  can  treat  any  arbitrary  boundary  condition.  Since 
the  finite  element  method  is  applicable  to  heat  transfer  analysis,  fluid  flow,  diffusion 
or  stress  analysis,  some  of  the  standard  structural  programs  can  be  used  in  the  analysis 
of  heat  transfer.  Also,  the  results  from  the  analysis  are  compatible  as  input  to  a 
finite  element  stress  analysis  program  (1).  When  properly  formulated,  the  finite  element 
method  gives  rise  to  a sparse  matrix  which  can  be  solved  very  efficiently  (PCCTSS-thermal 
program  uses  only  about  1 to  3 minutes  of  computer  time  on  the  Pitt  DEC  system).  Also, 
the  accuracy  of  solution  can  be  improved  by  increasing  the  number  of  elements  or  by  in- 
creasing the  number  of  nodes  per  element. 

In  this  paper,  the  application  of  the  finite  element  method  to  the  heat  transfer  analysis 
in  continuous  casting  is  outlined.  The  treatment  of  various  boundary  conditions  are 
explained.  The  results  from  the  analysis  are  compared  to  the  published  data  on  heat 
transfer  in  continuous  casting. 


Finite  Element  Analysis 

The  etepa  involved  in  the  finite  element  heat  transfer  analysis  are  the  following: 

1. )  Idealisation  of  the  casting 

2. )  derivation  of  the  element  stiffness  matrix  by  transforming  the 

governing  partial  differential  equation 

3. )  assembly  of  equations  and  application  of  proper  boundary  conditions 

4. )  solution  of  equations  for  nodal  temperatures 

5. )  computation  of  the  relevant  quantities  like  freezing  front  location, 

heat  flux,  etc. 

Each  of  these  steps  will  be  briefly  described  in  the  following  paragraphs. 

Idealisation  of  the  Casting 

The  cylindrical  ingot  is  divided  into  a number  of  triangular  ring  elements  as  shown  in 
Figure  1.  Two  of  these  triangular  ring  elements  grouped  together  give  the  quadrilateral 
ring  element  which  is  the  basic  element  used  in  the  present  analysis.  An  element  is  de- 
fined by  the  nodes  (e.g. , i,  J,  k in  Figure  lb)  and  nodal  temperatures.  The  casting  is 
considered  to  be  axisynnetric  and  hence  temperature  along  airy  of  the  circular  nodes  is 
constant. 

Derivation  of  the  Element  Stiffness  Matrix 


The  temperature  field  within  a triangular  element  may  be  expressed  in  terms  of  the  nodal 
temperatures  by  means  of  suitable  interpolation  functions. 

T « K1Ti  + NjTj  + 

* {Hj}  £T£}  CD 

where  (N^}  is  the  interpolation  function  and  (T^  is  the  column  matrix  of  nodal 
temperatures.  The  differential  equation  for  the  heat  flow  in  continuous  casting 
(axisymoatric)  is  given  by 

VkVT  -v  pCp^  = ?Cp  ||  (2) 

where  k is  the  thermal  conductivity,  v is  the  casting  rate,  p is  the  density,  Cp  is  the 
specific  heat,  and  9 is  the  time.  The  coordinate  system  is  regarded -as  fixed  with 
respect  to  the  caster  and  the  second  term  on  the  left  of  equation  (2)  accounts  for  heat 
transfer  by  mass  flow.  For  steady  state  conditions,  ^ = 0 and  equation  (2)  reduces  to 
Vk7T  - vpCp^|  » 0 (3) 

The  general  boundary  conditions  encountered  in  continuous  casting  are 

T * Tj  (temperature  specified)  (4) 

k^|  + q hT  * 0 (5) 

(heat  flux  or  convection) 

*|n  * - T*)  « 0 (6) 

(radiation) 

where  n is  the  normal  to  the  surface,  Tg  is  the  specified  surface  temperature,  q is  the 
heat  flux  from  the  boundary,  h is  the  heat  transfer  coefficient, r is  the  Stefan 


Boltzmann  constant,  £ is  the  emissivity  of  the  material,  and  T„  is  the  ambient 
temperature.  Applying  the  Galerkin  weighted  residual  method  to  equation  (3)  gives 
Vr  (<7k*T)  - vPcP$|)  dv  * 0 (7) 

where  dV  = 2TTrdrdZ 

Substituting  the  above  into  equation  (7)  and  integrating  by  parts  yields 

21prr  L + 

+ ypCpNi  drdZ  - \ [k^I  yfc 

♦ k^l  1ZNi]  dS  = 0 (*> 

where  lr  and  l£  are  the  direction  cosines  of  the  outward  normal  to  the  surface  and 

dS  is  the  differential  surface  element. 


In  assembling  these  for  the  entire  casting,  the  surface  integrals  will  cancel  on  neigh- 
boring elements  and  hence  will  give  a contribution  only  for  the  elements  which  have  an 
external  surface.  These  then  provide  a convenient  way  of  defining  the  boundary 
condition.  Equations  (4)  through  (6)  can  be  used  to  specify  the  boundary  contribution. 
For  example,  a convection  boundary  becomes 
Jc|I  + h(7-T.)  = 0 

Using  the  weighted  residual  method  this  transforms  to 
^ ^ dS  + h(T-Tj  N^dS]  * 0 

and  rewriting 

\ VW 

■ - ^ h(T-Tj  J^dS  (10) 

Equation  (10)  can  be  substituted  into  equation  (9)  for  the  element  having  an  external 
surface  subject  to  convection.  Heat  flux  or  radiation  boundary  conditions  can  be 
treated  in  a similar  manner. 


Substituting  aquation  (1)  to  equation  (9)  and  rearranging  terms  yields 

Chl  [Til  a W (11) 

where  [h]  is  the  element  thermal  stiffness  matrix,  is  the  nodal  temperatures,  and 

[q]  is  the  thermal  force  vector  (contribution  from  boundary  conditions). 

Assembly  and  Solution  of  Equations 

For  each  of  the  elements  in  the  body  a set  of  equations  similar  to  (II)  can  be  written. 

A number  of  elements  may  share  one  node  and  the  temperature  of  the  node  should  be  unique. 
This  condition  is  satisfied  by  assembling  all  the  equations  into  a global  equation. 
Assembly  yields  the  set  of  equations 

[H]  [9]  * [Q]  (12) 

where  (H]  is  the  global  thermal  stiffness  matrix,  (9}  is  the  vector  of  all  nodal 
temperatures  in  the  body  and  (Ql  is  the  thermal  load  vector  (contributions  from  all 


boundaries).  Since  both  [H]  and  [Q]  are  known,  equation  (12)  may  be  solved  to  get  the 
nodal  temperatures.  The  solutions  may  be  done  by  any  of  the  standard  methods  like 
matrix  inversion  methods,  Gaussian  elimination  method  or  iteration  methods  like  Gauss- 
Seidal  method.  Matrix  inversion  methods  consume  considerable  computer  time  when  a large 
number  of  equations  are  to  be  solved.  Gaussian  elimination  methods  give  an  exact  solu- 
tion except  for  roundoff  errors.  In  the  present  analysis,  the  Gauss-Seidal  iteration 
method  is  used.  This  avoids  the  problem  of  accumulation  of  roundoff  errors  and  converges 
quickly.  To  conserve  the  use  of  computer  .memory,  a sparse  matrix  method  is  used  storing 
only  the  nonzero  terms.  The  stiffness  matrix  obtained  here  is  nonsymmetric  (unlike  most 
of  the  structural  cases  where  the  matrix  is  symmetric)  and  hence  the  standard  solution 
schemes  for  solving  symmetric  matrices  cannot  be  used  here. 


Boundary  Conditions 


The  program  was  designed  so  that  it  will  take  any  boundary  condition.  However,  it  is 
of  interest  to  look  at  some  of  the  standard  boundary  conditions  encountered  in  continuous 
casting.  These  standard  boundary  conditions  were  used  in  getting  the  results  presented 
later. 

Center  of  Casting 

This  is  treated  as  an  insulating  boundary  condition 


'br 


= 0 


‘center 

This  is  the  natural  boundary  condition  in  finite  element  heat  transfer  analysis. 

Top  of  the  Casting 

To  "simulate"  a refractory  header  used  in  D.C.  casting,  the  top  ltver  of  elements  is 
fixed  at  the  pouring  temperature  [equation  (4)  with  ^pouring}* 

Bottom  Boundary 

The  bottom  boundary  is  considered  conducting.  The  total  heat  flux  is  by  conduction 
and  due  to  mass  movement. 

Circumferential  Boundary 

The  circumferential  boundary  consists  of  mold,  spray,  and  (in  the  case  of  steel) 
radiant  zones. 

Mold 

The  mold  is  considered  to  consist  of  a zone  of  good  contact  and  a zone  where  an  air 
gap  forms.  In  the  zone  of  good  contact  the  total  resistance  to  the  heat  flow  (.  -Tg00d ) is 
the  sum  of  the  resistances  due  to  a thin  film  of  oil  between  the  mold  and  the  casting 
(Uoii)  the  thickness  of  the  mold,  and  the  interface  between  the  mold  and  the  mold  cool- 
ing water:  x 

J m !1  J. 

good 


mold 


- 3 uoii + rr:  * t 


mold 


water 


- -i&msr  . 


where  is  the  thickness  of  the  mold,  k^^  is  the  conductivity  of  the  mold  material, 

and  hwater  is  the  heat  transfer  coefficient  at  the  mold/ water  interface.  As  the 

solidified  layer  increases  in  thictaess  and  cools  down,  it  develops  enough  strength  to 

withstand  the  metallostatic  head,  pulls  away  from  the  mold  due  to  the  solidification 

shrinkage,  and  an  air  gap  forms  between  the  mold  and  the  solidified  shell.  The  jpp  has  a 

strong  influence  on  the  total  heat  transferred  in  the  mold  which  in  turn  has  an  influence. 

on  gap  formation.  In  the  present  analysis,  the  complex  problem  is  treated  by  iteration 

(23).  As  will  be  seen  later,  this  procedure  gives  results  very  close  to  what  is  obtained 

experimentally.  The  thermal  resistance  from  the  shell  surface  to  the  mold  cooling  water 

in  the  presence  of  an  air  gap  is  given  by 

U i + U . 

gap  . good 

♦ <re(T4-T4) 

a*gap 

where  U_  is  the  total  thermal  resistance  from  the  skin  surface  to  the  mold  coolin'* 
gap 

water,  k^ap  is  the  thermal  conductivity  of  the  gas  in  the  gap,  agggp  is  the  thickness 
of  the  gap,  T is  the  temperature  of  the  surface  of  the  skin  and  Tm  is  the  temperature  of 
the  mold  surface.  Other  terms  are  defined  earlier. 

Spray  Zone 

In  the  spray  zone  the  heat  transfer  from  the  surface  of  casting  to  the  spray  water  can  be 
expressed  by  the  relation 

^ “ ^spray^T”Tspray^ 

where  T is  the  temperature  of  spray  and  h^^y  is  the  heat  transfer  coefficient  in 

the  spray,  hgp^y  can  be  obtained  from  experiments  such  as  those  reported  by  Mizikar 
(17).  Special  techniques,  such  as  delayed  quench,  may  be  simulated  by  varying  the  value 

^spray* 

Radiant  Cooling  Zone 

This  is  present  only  in  continuous  casting  of  steel.  The  heat  transfer  in  this  case  can 
be  given  by 

where  the  effective  heat  transfer  coefficient  in  this  zone,  h^,  is  given  by 
hrad  * <T£  (T2*  T2)  (T  + TJ 

Material  Properties 

The  program  has  the  capability  for  handling  fixed  and  variable  material  properties.  The 
variation  of  thermal  conductivity  and  density  with  temperature  have  been  extrapolated/ 
interpolated  from  known  values.  The  latent  heat  of  solidification  was  included  in  the 
heat  content  of  the  material  and  distributed  throughout  the  frequency  range  in  proportion 
to  the  solid  fraction  formed.  Both  equilibrium  and  nonequilibrium  solidification  modes 
can  be  handled  by  the  program.  The  convection  in  the  liquid  pool  was  handled  by  increas- 


ing  the  thermal  conductivity  of  the  liquid  metal  as  suggested  by  Mizikar  (17). 


Results 

In  order  to  test  the  accuracy  of  the  present  analysis,  the  results  obtained  from  the 
analysis  were  compared  with  those  reported  in  the  literature  for  laboratory  or  industrial 
tests.  Figures  2 and  3 show  the  comparison  of  the  model  predicted  temperature  profile 
with  those  reported  by  Adenis,  Coats  and  R*gone  (16)  for  magnesium  D.C.  casting.  These 
investigators  cast  thermocouples  on  the  surface,  center  and  2 inches  from  the  surface  of 

the  casting  and  recorded  the  temperature  profile  for  two  mold  diameters  (is. 5 and  9.5  cm.). 
The  predicted  profiles  are  found  to  be  very  close  to  those  obtained  experimentally. 

Table  1 compares  the  sump  depth  reported  by  Adenis,  Coats  and  Ragone  for  magnesium  alloy 
AZ-80A  D.C.  castings  with  those  obtained  from  the  present  analysis.  Again,  the  agreement 
is  found  to  be  very  good.  Since  the  magnesium  alloy  has  a very  high  thermal  conductivity 
and  is  cast  at  very  low  speed,  the  axial  heat  conduction  is  very  important  in  this  case. 
This  good  agreement  between  the  analysis  and  experimental  values  could  not  have  been 
obtained  if  the  axial  heat  conduction  were  ignored. 

Figure  4 shows  the  comparison  of  the  solidification  profiles  predicted  by  the  model  with 
those  reported  by  Weinberg,  et_al.(13),  from  tests  done  at  Western  Canada  Steel  Company. 
These  investigators  have  added  radioactive  tracers  to  the  pool  and  determined  the  solidi- 
fication profile  from  penetration  of  the  tracer  in  the  casting  cross  section.  The  agree- 
ment is  very  good  between  the  shell  thicknesses  predicted  from  the  analysis  and  the  pro- 
files obtained  from  industrial  tests.  Figure  5 and  Figure  6 compares  the  predicted  and 
measured  (13)  solidification  profile  for  plain  carbon  steel  and  stainless  steel  slab. 

In  measuring  these  profiles,  the  investigators  had  encountered  problems  in  penetration 
of  radioactive  gold  to  the  pool.  Thus,  the  profile  shown  as  measured  may  be  slightly 
in  error. 

Table  2 compares  the  temperature  of  the  skin  on  exit  from  the  mold  predicted  by  the  model 
and  those  obtained  experimentally  by  Gautier,  et  al.  (19).  These  investigators  have  re- 
corded the  temperature  of  the  casting  on  exit  from  the  mold  by  channelling  light  with 
light  pipes  from  the  billet  face  to  a two  color  pyrometer.  They  found  that  the  temper- 
ature of  the  four  faces  of  the  billet  fluctuated  a great  deal  during  the  experiment. 

The  values  shown  in  the  table  are  the  average  of  the  temperatures  recorded  at  the  four 
faces  of  the  billet.  The  analysis  is  found  to  predict  the  temperature  of  the  billet  on 
exit  from  the  mold  with  reasonable  accuracy. 


TABLE  1 


Comparison  of 

Predicted  and 

Measured  (16) 

Sump  Depth  of  Magnesium  Alloy  AZ80A  Billet 

3111«t  dia 
in 

Casting  Speed 
in/  min 

Mold 

length 

in 

Calculated 
Sump  depth 
in 

Measured 

Sump 

depth  in 

16 

2.1 

9.3 

18.17 

14.0 

16 

2.1 

16.3 

18.44 

18.3 

16 

3.0 

16.3 

21.60 

22.8 

TABLE  2 

Comparison  of  Predicted  and 

Measured  (19) 

Surface  Temperature  of  Plain  Carbon  Steel  Billet  on  Exit 

from  Mold 

Cast 

no 

Steel  Temp 
in  Distributor 
°C 

Withdrawal 
Speed 
cm/  sec 

Measured 
Billet  Surface 

**? 

Calculated 
Billet  Surface 
tSmp 

: 85  section 

42771 

1520 

4.9 

1237 

1231 

13400 

1503 

5.1 

1234 

1224 

13401 

1503 

5.9 

1228 

1241 

13402 

1S36 

5.1 

1210 

L287 

x 105 

section 

13761 

1505 

3.5 

1093 

1220 

1376  3 

1515 

3.3 

1089 

1220 

x 120 

section 

42803 

1508 

3.0 

1100 

1213 

13427 

1545 

2.7 

1106 

1241 

13428 

1546 

3.2 

1054 

1259 

Conclm  ions 


Tha  thermal  analysis  was  found  to  be  applicable  to  D.C.  casting  of  nonferrous  metals 
and  the  continuous  casting  of  steel.  The  comparison  of  the  results  obtained  from  the 
model  with  the  reported  experimental  values  indicate  that  the  agreement  is  very  good. 

The  analysis  provides  a powerful  tool  for  the  study  of  continuous  casting  process  when 
used  alone  or  along  with  a stress  analysis  model  described  in  another  paper. 
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FIG  I« 


Idealisation  of  the  casting. 


FIG  lb 

Basic  element  used  in  the  analysis 
i,  J and  k are  the  nodes  and  T £ , Tj 
and  are  the  nodal  temperatures. 
Notice  that  nodes  are  circuLar  here. 


Tempera  turu 


FIG  2 

Comparison  between  calculated 

( ) and  experimental  ( ) 

(16)  temperature  profiles  for 
magnesium  alloy  AZ80A  billet, 
billet  of  dia  16  in;  mold 
length  9.5  in;  casting  speed 
2.1  in/min. 


Comparison  between  calculated 

( ) and  experimental  ( ) 

(16)  temperature  profile  for 
magnesium  alloy  AX30A  billet, 
billet  dia  16  in;  mold  length 
16.5  in;  casting  speed  2.1 
in/min. 
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Comparison  between  calculated 

( ) and  measured  (•)  (13) 

solid iftcat ton  profile  14x14 
cm  low  carbon  steel  billet 
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FIG  S 

Comparison  of  the  solidifica- 
tion profile  obtained  from 
the  present  analysis  with 
those  reported  by  La it,  et 
tl.  (13)  for  13.3x13.3  cnT 
low  carbon  steel  billet 

analysis  3cir/sec; 

—analysis  Scm/sec. 


Shell  thickness  (cm) 


FIG  6 

Comparison  between  the  predicted 
(— ) and  measured  (•)  (13) 
solidification  profile  for  stain- 
less steel  slab. 
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SIMULATION  OF  THERMAL  STRESSES  IN 
CONTINUOUS  CASTING  USING  A FINITE  ELEMENT  METHOD 


J.  Mathew  and  H.  D.  Brody 
Metallurgical  and  Materials  Engineering 
University  of  Pittsburgh 


Abstract 

This  paper  presents  a three  dimensional  analysis  of  steady  state  thermal  stresses 
developed  in  the  continuous  casting  of  cylindrical  sections.  The  analysis  takes  into 
account  the  elastic,  plastic  and  steady  state  creep  deformations  of  the  casting.  It  can 
take  into  account  virtually  all  boundary  conditions  encountered  in  continuous  casting 
(friction  in  the  mold,  metallostatic  head,  etc.)  and  is  applicable  to  both  ferrous  and 
nonferrous  continuous  casting.  The  results  from  the  analysis  are  compared  to  the 
published  residual  stress  measurements  made  on  Aluminum-12%Si  alloy  (A132).  The  model 
is  very  useful  in  understanding  the  influence  of  the  operating  parameters  on  the  hot 
tearing  and  cold  cracking  of  continuous  castings. 

Introduction 

One  of  the  continuing  problems  in  the  continuous  casting  of  ferrous  and  nonferrous  alloys 
is  the  incidence  of  cracking  in  the  cast  section.  For  example,  aluminum  alloys  (2024  and 
7075)  are  cast  at  extremely  low  speeds  to  avoid  cracking.  Higher  casting  rates  would 
be  desirable  for  reasons  of  productivity,  surface  finish,  etc.  The  main  cause  for  the 
formation  of  cracks  is  the  development  of  thermal  stresses  in  the  casting  during  solidi- 
fication and  subsequent  cooling.  Hence,  to  understand  the  influence  of  operating  para- 
meters on  the  sensitivity  of  the  casting  to  cracking,  it  is  helpful  to  compute  the 
thermal  stresses  developed  in  the  casting.  This  paper  describes  a simulation  program 
developed  at  the  University  of  Pittsburgh  to  compute  the  thermal  stresses  in  continuous 
casting.  When  used  along  with  a heat  transfer  model  of  continuous  casting  (1),  the 
analysis  provides  a powerful  tool  for  understanding  the  origin  of  various  thermal 
stress  related  defects  in  continuous  casting. 

The  analysis  presented  in  this  paper  is  applicable  to  ferrous  and  nonferrous  continuous 
castings  of  cylindrical  section.  It  is  a three  dimensional  axisymmetric  analysis  which 
takas  into  account  the  elastic,  thermal,  plastic  and  creep  deformations  of  the  material. 
The  variation  of  mechanical  properties  of  the  material  (elastic  modulus,  Poisson's  ratio, 
coefficient  of  linear  expansion,  yield  stress,  plastic  modulus  and  creep  rate  with 
temperature)  can  be  taken  into  account  in  this  analysis.  The  effect  of  body  forces. 


friction  forces  in  the  mold,  metallostatic  heed,  pinch  roll  pressure,  and  withdrawal 
forces  can  be  considered.  An  automatic  feature  of  the  program  is  to  set  the  thermal 
boundary  condition  in  mold  based  on  computing  the  position  of  formation  and  the  width 
of  the  air  gap.  As  the  computation  of  stresses  in  the  casting  can  be  continued  until  the 
casting  attains  room  temperature,  the  present  analysis  can  be  used  to  compute  the 
residual  stresses  in  the  casting.  The  residual  stresses  can  be  measured  experimentally, 
by  techniques  such  as  Sach's  boring  method.  The  simulation  model  can  be  used  to  predict 
the  state  of  stress  during  solidification  (which  is  very  difficult  to  measure)  after 
calibration  with  measurements  on  the  casting  of  room  temperature  residual  stresses. 

In  this  paper,  we  shall  briefly  discuss  the  finite  element  analysis  of  thermal  stresses 
in  the  casting.  The  results  from  the  analysis  are  then  compared  to  the  residual  stress 
measurements  made  on  Aluminum  A132  alloy.  The  utility  of  the  simulation  model  to  study 
the  hot  tearing  and  cold  cracking  of  the  continuously  cast  sections  will  be  indicated. 

Analysis 

The  idealisation  of  the  casting  is  similar  to  the  one  used  for  the  heat  transfer  analysis 
(1).  The  basic  element  is  a triangular  ring  element  with  three  nodes.  The  displacement 
field  (u)  within  this  element  is  expressed  in  terms  of  nodal  displacements  (5)  by  means 
of  suitable  interpolation  functions  [N]  (2): 

(u)  * [N]  [a]  (1) 

From  the  strain-displacement  relations,  the  strain  within  the  element  can  be  computed  as 

c«  ■ :Bi  {«}  <2> 

where  (£}  is  the  column  matrix  of  element  strains  and  [B]  is  the  strain  displacement 
matrix.  From  the  stress-strain  relations  for  the  material  we  get  the  relation 

frr}  * [D]  ([»}  - [%))  + & (») 

where  (3*}  is  the  vector  of  element  stresses.  [D]  is  the  stress-strain  matrix  which 
depends  on  the  properties  of  the  material  and  the  type  of  deformations  considered.  For 
the  axisymmetric  problem  considered  here,  the  strain  and  stress  vectors  are: 

tel  M 

*(t}  * w 7 awl  7 (4) 

1®tzJ  \^z  J 

For  instance,  if  elastic  analysis  is  carried  out,  the  matrix  [D]  can  be  derived  from 
generalised  Hooke's  law.  If  plasticity  or  creep  is  to  be  considered,  the  analysis 
becomes  nonlinear  as  the  matrix  [D]  now  depends  on  the  strain  [g] . The  technique  of 
handling  plasticity  and  creep  analysis  will  be  discussed  later.  In  equation  (3),  (£^3 
is  the  initial  strain  matrix  (thermal  deformation  is  a very  good  example  of  the  presence 
of  an  initial  strain)  and  is  the  initial  stress  matrix.  The  equivalent  nodal  forces 

that  are  to  be  applied  at  the  nodes  can  be  obtained  by  the  application  of  the  virtual 
work  theorem. 


W&-:  ■ 


(5) 


to*  * co*  ci!  * tnJo  * tr£o  * mj 

[H*  are  the  element  nodal  forces,  [k]*  is  ths  alamant  stiffnass  matrix,  [F]*  , [F^q, 
and  [FI*  are  tha  equivalent  nodal  forces  due  to  initial  strains,  initial  stresses,  and 
body  fortes,  respectively. 

Tba  equilibrium  relations  of  tha  stress  analysis  are  invoked  in  the  derivation  of  the 
element  stiffness  matrix.  The  compatibility  conditions  are  enforced  by  constraining  the 
displacements  at  the  nodes  of  the  element  to  be  unique.  This  is  achieved  by  assembling 
tha  element  stiffness  matrices  into  a global  stiffness  matrix  for  tha  casting  and  defin- 
ing a global  matrix  of  nodal  displacements.  This  gives 

[K]  [5)  = [F]  _ (6) 

Hera  [K]  is  the  global  stiffness  matrix,  [u]  is  the  global  displacement  matrix  and  [F]  is 
the  force  matrix  for  the  casting.  Equation  (6)  can  be  solved  for  tha  casting  nodal 
displacements.  The  solution  may  be  carried  out  by  any  of  the  standard  methods  like 
elimination  methods  or  iteration  methods.  (Matrix  inversion  methods  are  not  recommended 
as  they  involve  considerable  computer  time.)  In  the  present  analysis,  the  Gaussian 
elimination  method  was  used  for  the  solution  of  the  equations. 

From  the  casting  nodal  displacements,  the  strains  in  the  individual  elements  are  computed 
using  equation  (2).  The  stresses  in  the  element  are  computed  from  the  strains  using 
equation  (3).  From  this,  the  principal  stresses  and  strains  are  computed. 

Handling  of  Plasticity,  Creep  and  Thermal  Strains 

The  handling  of  plastic,  creep  and  thermal  deformations  in  .he  analysis  needs  special 
attention.  The  total  strain  can  be  considered  to  consist  of  elastic,  plastic,  creep 
and  thermal  components. 

^total  * ^elastic  + ^thermal  + Aplastic  * ^creef  ^ 

The  elastic  strains  are  related  to  the  stresses  by  the  Generalised  Hooke's  Law, 

Wi.  * **  <•> 

where  E is  the  modulus  of  elasticity  of  the  material.  The  thermal  deformations  of  the 
material  affects  only  the  dilatation  and  will  not  influence  the  shear  strains.  The 
Initial  strain  dua  to  tha  therrtml  strain  thus  become 


where  a is  the  linear  coefficient  of  expansion  of  the  material  and  aT  is  the  change  of 
temperature  from  the  stress  free  state. 

The  plastic  strains  are  related  only  to  the  deviatoric  stresses  and  they  are  not 
influenced  by  the  hydrostatic  stress.  [Plastic  strains  cannot  result  in  any  volume 


change.]  The  relation  between  the  incremental  plastic  strains  and  the  deviatoric 


' be  related  by  Prandtl  Reuse  relations 

(2<rr  - <re  - c rz) 

(2<r©  ’ <re  - <r7i) 

(2<rz  ‘ <Tr  ‘ aV 

AcP  «-L 

A*rZ  E„ 
P 

In  equation  (10),  E is  defined  as 
E -IS 


where  q-  is  the  equivalent  stress  defined,  by 

e *zzrV,(o-r-cre)2  + (cr0-<rz)2  + (c rr-<rz)2  + *4z 
V 2 


(10) 


(ii) 


and  Is  the  incremental  plastic  strain.  For  the  elements  which  have  stresses  above 

the  yield  stress  of  the  material  (determined  by  Mises'  yield  criterion),  the  matrix  [D] 
in.  equation  (3)  is  modified  into  a plastic  stress-strain  matrix  [Dp].  In  order  to  con- 
serve computer  time,  the  element  stiffness  matrix  is  left  unaltered.  This  can  be 
achieved  by  writing 

[Dp]  = [D]  + [O'] 

where  [Dp]  is  the  actual  plastic  stress-strain  matrix,  £D]  is  the  elastic  stress-strain 
matrix  and  [D']  is  the  difference  between  [Dp]  and  [D].  The  relation  between  stress  and 
strain  can  be  written  as 


w * rDPi  w 

* [D]  (£}  * [D']  ce} 

CCt]  + C£H  d2) 

Thus,  the  plastic  strains  result  in  initial  strain  . which  do  not  alter  the 
stiffness  matrix  of  the  element.  Only  the  force  vector  is  affected.  The  stiffness 
matrix  may  be  triangularised  after  the  assembly  process.  The  subsequent  solutions 
will  affect  only  the  back  substitution  process.  This  method  saves  considerable  amount 
of  computer  time  as  in  the  Gaussian  elimination  method  triangularisation  of  the  matrix 
consumes  most  of  the  computer  time. 


The  creep  strains  are  handled  in  an  identical  manner.  The  creep  rate  is  related  to  the 
equivalent  stress  by  the  creep  law  of  the  material.  From  the  creep  strain  and  equivalent 
stresses,  the  individual  components  of  strain  are  computed  using  equation(lO).  The 
creep  strains  are  handled  as  initial  strains  and  the  equivalent  nodal  forces  are 


computed.  The  new  casting  nodal  displacements  are  computed  by  back  substitution. 

Boundary  Conditions 

The  simulation  program  was  designed  to  handle  the  widest  range  of  boundary  conditions. 

The  boundary  conditions  used  to  obtain  the  results  presented  later  are  discussed  in  this 
section. 

Bottom  of  the  Casting 

The  boundary  is  considered  to  be  fixed  axially;  but  free  to  move  radially. 

Center  of  Casting 

The  casting  is  fixed  radially,  but  is  free  to  move  axially. 

Top  of  the  Casting 

Above  a certain  temperature  a liquid  plus  solid  mixture  behaves  like  a liquid  and  below 
that  temperature  the  mixture  behaves  like  a solid.  This  temperature  is  termed  a coherency 
temperature.  To  account  for  the  mushy  rone  in  the  casting,  the  elements  which  are  above 
the  coherency  temperature  are  removed  and  on  the  surface  of  the  resulting  boundary  the 
oetallostatic  head  is  applied. 

Circumferential  Boundary 

This  consists  of  mold  zone  and  the  spray.  In  the  spray  zone,  the  boundary  is  free.  In 
the  mold  the  elements  may  or  may  not  be  constrained  radially  depending  on  whether  an  air 
gap  is  present.  In  the  program,  the  mold  boundary  condition  may  be  specified  or  it  may 
be  generated  by  the  program  itself.  In  the  latter  case  the  casting  surface  is  allowed  to 
move  inward  or  outward  except  that  the  mold  diameter  is  considered  a limiting  dimension. 

Results 

Before  discussing  the  application  of  the  analysis  to  continuous  casting,  two  methods  of 
checking  the  accuracy  of  the  analysis  are  discussed.  First,  the  accuracy  of  the  displace- 
ment field  was  checked  by  applying  the  analysis  to  certain  ideal  cases  where  closed  form 
analytical  solutions  are  available  in  the  literature.  Both  thermal  stress  and  creep 
analyses  were  checked.  The  accuracy  of  the  model  when  applied  to  continuous  casting  was 
checked  by  comparing  the  model  predicted  stress  values  to  published  residual  stress 
measurements  made  on  continuous  casting. 

The  thermal  stress  analysis  model  (elastic)  was  checked  by  applying  the  analysis  of  a 
hollow  thick  cylinder  with  a steady  state  temperature  distribution.  The  elastic  thermal 
stresses  developed  in  the  cylinder  were  computed  and  compared  to  the  results  of  an 
analytical  solution  to  this  problem  given  by  Johnson  and  Mellor  (3).  Figure  1 compares 
the  model  predicted  stress  values  with  the  exact  solution  for  a hollow  cylinder  of  inside 
and  outside  radii  of  10  and  20  inches,  respectively,  with  the  inner  surface  maintained  at 
100°C  and  outer  surface  maintained  at  200°C.  The  agreement  is  found  to  be  very  good 
between  the  analytical  results  and  the  exact  solution. 


The  steady  state  creep  analysis  was  checked  by  applying  the  model  to  the  case  of  a 
hollow  thick  cylinder  with  a uniform  pressure  inside  the  cylinder.  The  elastic  and 
steady  state  creep  stresses  have  been  calculated  by  Odquist  (4)  for  this  case  when  the 

steady  state  creep  rate  is  given  by  the  law 

1 -18  4.4 

6 = 6.4  x 10  ^ 

where  £ is  the  creep  rate  and  a*a  is  the  equivalent  stress.  The  model  was  applied  to 
the  case  of  a cylinder  of  inner  radius  0.16”  and  outer  radius  0.25”  with  a uniform 
inner  pressure  of  365  psi.  The  elastic  and  creep  stresses  are  shown  in  Figure  2.  Again, 
note  that  the  agreement  between  the  exact  solution  and  the  model  predicted  values  is  very 
good. 


Using  Sachs  boring  technique,  Roth,  et  al.  (5)  measured  the  residual  stresses  in  6.5  inch 
and  5.7  inch  diameter  Al-12%Si  alloy  ingots.  The  ingots  were  cast  at  speeds  of  8.95  and 
17.1  in/min.  The  6.5  inch  diameter  ingot  was  found  to  have  a residual  longitudinal 
tensile  stress  of  26,000  psi  at  the  axis  when  cast  at  these  rates.  This  was  found  to 
decrease  toward  the  surface  and  changed  to  compressive  stresses  1.14  inches  below  the 
surface,  the  stresses  was  -14,000  to  -18,000  psi.  Roth  concluded  that  for  the  two  ingot 
sizes  tested,  the  rate  of  casting  had  no  marked  effect  on  the  room  temperature  residual 
stresses.  But  the  cracking  tendency  was  found  to  increase  with  the  increase  in  casting 
speed. 

The  present  analysis  was  used  to  simulate  the  casting  of  A132  alloy  ingots  6.5  inches 
in  diameter.  The  axial  stresses  obtained  from  the  present  analysis  are  shown  in  Table  1-2 
for  casting  speeds  of  3.95  iry'min  and  17.1  in/min.  The  results  are  summarized  in  Table  3 
along  with  the  experimental  results  obtained  by  Roth,  et  al.  Notice  that  the  agreement 
is  fairly  good  between  the  model  predicted  residual  stresses  and  those  obtained  experi- 
mentally. The  present  analysis  also  predicts  quite  accurately  the  reversal  of  the 
stresses  from  tensile  to  compressive  near  the  surface. 

The  cracking  generally  originates  near  the  end  of  solidification  and  hence  the  stresses 
during  the  solidification  of  the  alloy  are  of  interest.  To  date  these  have  not  been 
obtained  experimentally.  Analyses  of  type  presented  here  are  necessary  to  estimate 
stresses  in  the  casting  during  solidification. 

Since  the  temperature  changes  from  the  top  to  bottom  of  the  casting,  the  absolute  stresses 
as  given  in  Tables  1 and  2 do  not  give  much  insight  as  to  the  severity  of  the  stress 
level.  Thus,  the  stress  values  were  normalised  by  dividing  by  the  yield  stress  of  the 
material  at  the  temperature  of  the  element.  The  maximum  normalised  stress  (maximum  of 
hoop,  radial  or  axial  stress)  at  different  levels  of  continuously  cast  A132  is  given 
in  Tables  4 and  5.  (The  temperatures  at  the  points  where  the  stress  is  computed  are 
tabulated  below  the  stress  values.)  Notice  thct  for  the  casting  conditions  simulated 
here,  the  maximum  normalised  stress  appears  just  at  the  solidification  front.  The  value 
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Computed  Axial  Stresses  in  a Cylindrical  Al-12%Si  Alloy  (A132)  Casting  of  Diameter 

6.S  inches  Cast  at  8.95  in/min 
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TABLE  2 

Computed  Axial  Stresses  in  a Cylindrical  Al-12%Si  Alloy  (A132)  Casting  of  Diameter 

6.5  inches  Cast  at  17.1  in/min 
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Computed  Maximum  Normalised  Stress  and  Temperature  Distribution 
in  an  A132  ALloy,  6.5  in.  dia.  Casting  Poured  at  8.95  in/min 
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Computed  Maximum  Normalised  Stress  and  Temperature  Distribution 
in  an  A132  Alloy,  6.5  in.  dia.  Cisting  Poured  at  17.1  in/min 
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Is  2.15  for  the  ingot  cast  at  8.95  in/min  and  4.97  for  the  ingot  cast  at  17.1  in/min. 
This  indicates  that  the  ingot  cast  at  17.1  in/min  is  snore  prone  to  hot  tearing  than  the 
one  cast  at  8.95  in/min  (as  observed  by  Roth).  Since  the  maximum  stress  appears  at 
the  axis  of  the  casting,  the  crack,  if  formed,  will  be  axially  symmetric  (ghost  or  star 
cracks). 

Summary 

The  analysis  presented  in  this  paper  is  very  useful  in  studying  the  effect  of  operating 
parameters  on  continuous  casting.  As  the  stresses  in  the  solidification  range  as  well 
as  down  to  room  temperature  can  be  obtained,  it  can  be  used  to  simulate  both  hot  tearing 
and  cold  cracking  of  the  casting.  The  analysis  gives  the  principal  stresses,  maximum 
normalised  stresses  at  different  points  in  the  casting,  plastic  strain  in  the  element 
and  final  diameter  of  the  casting  after  cooling.  This  information  is  sufficient  to 
arrive  at  a criterion  for  hot  or  cold  cracking  of  the  casting.  It  is  applicable  to  any 
alloy  (aluminum  O.C.  casting  as  well  as  continuous  casting  of  steel).  The  accuracy 
of  the  analysis  can  be  checked  by  room  temperature  residual  stress  measurements.  Mold 
length,  heat  transfer  in  the  secondary  cooling  zone,  superheat,  etc.  are  very  important 
parameters  which  Influence  the  hot  tearing  tendency  of  the  material. 

The  agreement  between  the  hot  tearing  tendency  of  castings  as  predicted  from  the  present 
model  apd  the  practical  experience  in  the  D.C.  casting  of  Aluminum  6063,  5084  and  2014 
alloys  are  found  to  be  very  good.  The  model  is  being  used  at  present  for  a parametric 
study  of  the  hot  tearing  in  continuously  cast  sections. 
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^ This  report  sumoarlzes  work  on  composite  and  alloy  solidification 


presently  conducted  under  the  ONR  contract  and  includes  two  papers  on  the 
simulation  of  heat  flow  and  thermal  stresses  in  continuous  casting  presented 
at  the  "Conference  on  Computer  Simulation  for  Materials  Applications"  in 
Gaithersburg,  Md.  in  April  1976. 

A final  element  analysis  of  heat  transfer  in  continuous  casting  is 
presented  in  the  first  paper.  The  steady  state  conditions  of  the  solidifies 
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tion  of  a cylindrical  alloy  ingot  are  simulated.  The  analysis  takes  into 
account  both  radial  and  axial  heat  conduction  and  the  variation  of  alloy 
properties  with  temperature.  The  analysis  has  the  versatility  to  simulate 
any  casting  configuration  as  the  program  has  the  capability  to  treat  any 
arbitrary  boundary  conditions.  In  the  mold  region,  where  an  air  gap 
generally  forms,  the  position  of  gap  formation,  the  size  of  the  gap  and  the 
heat  transfer  through  the  ^p  by  conduction  and  radiation  are  computed. 

Values  predicted  by  the  model  are  in  close  agreement  with  published 
laboratory  and  industrial  data . _ 

k: 

The  second  paper  presents  a three  dimensional  analysis  of  steady  state 
thermal  stresses  developed  in  the  continuous  casting  of  cylindrical  sections. 
The  analysis  takes  into  account  the  elastic,  plastic  and  steady  state  creep 
deformations  of  the  casting.  It  can  take  into  account  virtually  all  boundary 
conditions  encountered  in  continuous  casting  (friction  in  the  mold,  metallo- 
static  head,  etc.)  and  is  applicable  to  both  ferrous  and  nonferrous  continuous 
casting.  The  results  from  the  analysis  are  compared  to  the  published 
residual  stress  measurements  made  on  aluminum- 12%  silicon  alloy  (A132 ) . 

The  model  is  very  useful  in  understanding  the  influence  of  the  operating 
parameters  on  the  hot  tearing  and  cold  cracking  of  continuous  castings. 
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